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Abstract 

We study and compare three coded schemes for single-server wireless broadcast of multiple description coded 
content to heterogeneous users. The users (sink nodes) demand different number of descriptions over links 
with different packet loss rates. The three coded schemes are based on the LT codes, growth codes, and 
randomized chunked codes. The schemes are compared on the basis of the total number of transmissions 
required to deliver the demands of all users, which we refer to as the server (source) delivery time. We 
design the degree distributions of LT codes by solving suitably defined linear optimization problems, and 
numerically characterize the achievable delivery time for different coding schemes. We find that including a 
systematic phase (uncoded transmission) is significantly beneficial for scenarios with low demands, and that 
coding is necessary for efficiently delivering high demands. Different demand and error rate scenarios may 
require very different coding schemes. Growth codes and chunked codes do not perform as well as optimized 
LT codes in the heterogeneous communication scenario. 
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1. Introduction 

1.1. Motivation 

In the past decade, the development of wireless networks has provided a fertile soil for popularization of 
portable digital devices, and the digital distribution of bulk media contents. This, in return, is stimulating 
new leaps in wireless communication technology. Today, devices retrieving digital video content in the air 
vary from 1080p HDTV sets to smartphones with 480x320- pixel screen resolution. Under assorted restraints 
in hardware, power, location, and mobility, these devices experience diverse link quality, differ in computing 
capability needed to retrieve information from received data streams, and request information of varied 
granularity. Consequently, a transmission scheme designed for one type of users may not be as suitable for 
another even if they demand the same content. 

Today's technology implements a straightforward solution of separate transmissions to individual users, 
that is, multiple unicasts. Nonetheless, the key question is whether it is possible to deliver all users' demands 
with fewer data streams and less traffic. Especially in transmitting bulk data through wireless channels or 
over other shared media, reducing the amount of traffic is vital for reducing collision/interference, which in 
turn will also positively affect the quality of the channel in use. An additional concern in the environment 
conscious world is, of course, energy. 

There is no surprise then that the problem of delivering more efficient service to a heterogeneous user 
community has attracted a great deal of both technical and academic interest. On the source coding side, 
layered coding (e.g.[l|, Q) and multiple description coding Q have been widely studied as solutions to 
providing rate scalability. In particular, with multiple description coding, a user is able to reconstruct a 
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lower-quality version of the content upon receiving any one of the descriptions, and is able to improve the 
reconstruction quality upon receiving any additional description. Thus, the users' content reconstruction 
Quality is commensurate with the quality of their connection. On the channel coding side, rateless codes 
4, 5], or fountain codes, can generate a potentially infinite stream of coded symbols that can be optimized 
simultaneously for different channel erasure rates as long as the users have uniform demands. 

In this work, we explore an achievable efficiency of serving users having heterogeneous demands while 
using a single broadcast stream. Whearas the information theoretical aspect of the problem is of interest 
and under investigation (see for example @ and references therein), we focus on three practical coding 
schemes and explore their suitability for the described communications scenario. Two important features 
that make codes suitable for such scenarios are (1) the ability to support partial data recovery and (2) the 
ability to efficiently adapt to different channel conditions. Based on these desirable features, we chose to 
investigate three candidates: LtQ, growthQ, and chunked d, Q codes. In this paper, we are particularly 
interested in the total number of source transmissions needed to deliver the demands of all users. This 
quantity determines the amount of required communication resources, and also translates to the amount 
of time required for delivery. In the streaming of temporally-segmented multimedia content that is delay- 
constrained, it is important for the source to finish transmitting a segment as soon as possible so as to 
proceed to the next one. Some performance measures of interest are addressed in (Io| . 

1.2. Main Results 

We compare three coding schemes and their variations for broadcasting to heterogeneous users: users 
suffer different packet loss rates and demand different amounts of data. The coding schemes discussed 
include: 

1. optimized LT codes, with or without a systematic phase, that is, one round of transmission of the 
original uncoded packets; 

2. growth codes; and 

3. chunked codes. 

We also compare these schemes to a reference scheme for the heterogeneous scenario based on time-shared 



broadcast of degraded message sets[ll(. We find that including a systematic phase is significantly beneficial 
towards delivering lower demands, but that coding is necessary for delivering higher demands. Different user 
demographics result in the suitability of very different coding schemes. Growth codes and chunked codes are 
not as suitable to this communication scenario as are optimized LT codes. 



1.3. Organization 

The rest of the paper is organized as follows. Section [2] introduces the model of wireless broadcasting 
of multiple description coded content to heterogeneous users. Section [3] introduces the coding schemes of 
interest, and provide theoretical characterization of code performance. In particular, the LT codes are 
optimized both with and without a systematic phase. In Section [4] we provide numerical and simulation 
results of the achievable server delivery time, and discuss the suitability of these coded schemes for broadcast 
to heterogeneous users. The last section concludes. 



2. System Model 

Consider a wireless single-hop broadcast network consisting of a source (server) node holding content for 
distribution, and I sink (user) nodes waiting to retrieve the content from the broadcast stream aired by the 
source, as shown in Figure [T] Suppose the content held at the source node is multiple description coded into 
N descriptions, using, for example, one of the coding schemes described in [l2[ or in [Hj], and each description 
is packaged into one packet for transmission. Each packet is represented as a binary vector of length b, and 
denoted with for j = 1, 2, . . . , N. A low-quality version of the content can be reconstructed once a user is 
able to recover any description. The reconstruction quality improves progressively by recovering additional 
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Figure 1: Network model 



descriptions, and depends solely on the number of descriptions recovered, irrelevant of the particular recovered 
collection. 

The source broadcasts one packet per unit time to all the sink nodes in the system. However, the sinks 
are connected to the source by lossy links, and at every sink, a packet either arrives at the sink intact or is 
entirely lost. Such an assumption is practical if we only consider data streams at the network level or higher. 
Under the multiple description coding assumption, a sink node can choose to demand a smaller number of 
descriptions rather than wait to collect all N. This may not only shorten its own waiting time but also 
reduce the system burden. The demand of a sink node is characterized by the number of descriptions it 
needs to reconstruct the content within the desired distortion constraint. 

Sink nodes are characterized by parameter pairs (zi, i = 1, 2, . . . , Z, describing their demands and link 
qualities: 

• Zi € {1/N,2/N, . . . , 1} is the fraction of content demanded by user i, that is, each node demands ZiN 
distinct descriptions. 

• a is the packet erasure (loss) rate on each link from the source to user i. A packet transmitted by the 
source fails to reach user i with probability a. All links are assumed to be memory less. 

We further assume no feedback from the sinks to the source except for the initial requests to register the 
demands and the final confirmation of demand fulfillment. We also assume that the packet erasure rates 
£iS on the links are known to the server. In wireless broadcast, feedback from multiple nodes increases the 
chance of collision and compromises the point-to-point channel quality. Hence, it is desirable to keep the 
feedback at a minimum. 

Definition 1. (General Delivery Time) The delivery time Tj of a sink is a random variable defined as the 
number of packets transmitted from the server until the fulfillment of the demand Zi of user i. The server 
delivery time T = max' =1 Ti is the number of packets transmitted from the server required to fulfill the 
demands of all the users. 

Since all users retrieve information from the same broadcast stream through channels of random erasures, 
Tj's and T are dependent random variables. Nevertheless, even their marginal probability distributions are 
not easy to characterize. Instead, in this paper, we normalize the delivery time by the total number of 
content packets N, and restrict our attention to either the expectation ti — E[Ti]/N or the asymptotics 
U = linijv->oo Ti/N, and study max' =1 ti as the (normalized) server delivery time. When ti = E[Ti]/N, 
max' =1 ti is a lower bound for E[T]/N, the normalized expected delivery time. For brevity we further abuse 
our notation and terminology and specialize the definition of delivery time to these lower bounds. The 
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specialization of a definition will be stated where appropriate. In the rest of this paper, we describe and 
analyze three randomized coding schemes for broadcast in a heterogeneous setting, optimize if appropriate, 
compare their performance, and discuss their suitability to our end. 

Note that with multiple unicasts, the total number of packets transmitted would be (in the mean) at 
least X^i=i T^ - ' W e will show in our results that it is possible to do better using broadcast. 



3. Three Coding Schemes for Broadcast to Heterogeneous Users 

In this section, we describe three coding schemes and characterize their delivery time performance for 
broadcast to heterogeneous users. All the coding schemes studied in this paper are applied at the packet 
level. Each coded packet is obtained as a random linear combination of the original packets, but either the 
number or the range of the packets involved in each linear combination is restricted. 



3.1. LT Coded Broadcast 

LT codes were first proposed by Michael Luby in 0] as a class of codes designed to overcome variations 
in channel quality via ratelessness. Unlike their rateless predecessors, the ordinary random linear codes, 
LT codes can be decoded by simple suboptimal decoders with little loss in performance. In [13], Sujay 
Sanghavi pointed out the inadequacy of the original LT codes for partial data recovery, and demonstrated 
the possibility of redesigning the LT codes. His focus was on scenarios with users having uniform demands. 
Here we extend the code optimization framework of 14[ to adapt the code to the heterogeneous scenario in 
which both the channel conditions and the demands vary among users. 



3.1.1. Overview of LT Encoding and Decoding 

The LT codes, as defined in Luby's original paper are binary. An LT encoder outputs a potentially 
infinite stream of random linear combinations of the N original content packets, and decodes with a belief- 
propagation decoder. Each coded packet is generated independently as follows: first, pick a degree d according 
to a distribution defined by its moment generating function 

P {N) (x)= Pl x+p 2 x 2 + ...+p N x N } (1) 

where pj = Prob[d = j]; then, select with equal probability one of the ( ,) possible combinations of d distinct 
packets from the N content packets, and form a coded packet of degree d by linearly combining the chosen 
packets over GF(2). The linear combination is represented by a coding vector c = [c\, c 2l . . . , c/v] T , a binary 
column vector of length N with exactly d non-zero components corresponding to the chosen packets. The 
generated coded packet is given by Y]j—x Cj£j (recall £jS are the binary vector representations of the content 
packets as defined in Section[2]) and it is of the same size b as the original content packets. The iV-bit coding 
vector c can be embedded into the coded packet and transmitted to the user. In wireless communications 
where data packets are generally of at most a few kilobytes' size, if N ~1K, it is costly to add an extra 
N bits. Instead, the source can embed into the packet the degree of the coding vector and the seed of a 
pseudo-random number generator that is in shared knowledge of both the source and all the sink nodes, to 
allow the sink node to regenerate the coding vector locally. 

Decoding is done by a belief propagation decoder [J, [5j. The decoder maintains a set called the ripple. 
The initial ripple is composed of all received coded packets of degree-1. In each decoding step, one coded 
packet in the ripple is processed by removing it from the ripple, substituting it into all the coded packets 
it participates in, which reduces the degrees of these coded packets. The coded packets whose degrees are 
reduced to 1 are added to the ripple, and the decoder continue to process another coded packet in the ripple. 
The decoding process halts when the ripple becomes empty. 
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3.1.2. Characterization of the Recoverable Fraction 

Each user i needs to recover z%N content packets from the received coded packets. 

Definition 2. (Recoverable Fraction) A fraction z is said to be recoverable by the belief propagation decoder 
if the ripple size stays positive (at least) until zN packets get decoded. 

Based on Theorem 2 in [15(, Corollary [2] gives the expected ripple size as a function of the recovered 
fraction of the content packets. We restate the part of Theorem 2 in [l5{ that concerns the expected size of 
the ripple as Theorem [1] 

Assume w ■ N coded packets have been collected and fed into an LT decoder, for some positive constant 
w. Let u ■ N be the number of decoded content packets, for a constant u 6 [0, 1). Let r^ N > (u) be the expected 
size of the ripple, normalized by N. 

Theorem 1. (Maatouk and Shokrollahi FZ1 . Thm. 2]) If an LT code of N content packets has degree distri- 
bution specified by the moment generating function P^ N \x) (see (GPj, then 

r W(u) = wu( f(% -u) + - In u) + o(— V (2) 
V w J \N J 

where p( N >'(x) stands for the first derivative of (x) with respect to x. 

We further extend Theorem [1] to the case where the number of collected coded packets is random to 
accommodate random packets losses over the channel. Assume the number of collected coded packets is now 
W ■ N, where W is a random variable with mean v. Denote the normalized expected ripple size as N — > oo 
as rw{u). Assuming that P( N '(x) converges to P(x) = X^iP* 2 " 1 as N — > oo, then the following holds: 

Corollary 2. 

r w (u) = u(vP'(l-u) + lnu). (3) 

PROOF. Although Theorem Q] (Theorem 2 in |15]) is stated for the case where the number of coded packets 
collected by the sink node is more than the total number of content blocks, i.e., w > 1, its proof suggests 
that the theorem also holds for any constant iu < 1. 
Take N — > oo on both sides of [2] we have 

r(u)= lim r {N) (u) = u(wP'{\ - u) +\nu). (4) 

Replace w in (Q} with W. Due to the linearity of the expected ripple size in W for given u and P, 

r w (u) = E u(wP'(l-u)+lnuj = u(vP'(l - u) + huj 

If we use the expected value as a rough estimate of the ripple size during the decoding process, we should 
have rw(u) > for u £ (1 — z, 1], Applying ([3]), we have 

vP'(l -u)+ In it > 0, Vue(l-z,l], (5) 

We next use (J5]) as a constraint to formulate an optimization problem for LT code degree distribution design 
to minimize the number of transmissions required to meet all sink demands. 

3.1.3. Minimizing Server Delivery Time by Degree Distribution Design 

We are concerned with the delivery time of LT coded broadcast in the asymptotic regime when N — > oo. 



Definition 3. (LT Delivery Time) For the LT coded scheme, the (normalized) delivery time ti is defined as 
the ratio of the number of transmissions required to fulfill the demand Zi of user i to the the total number 
N of content packets, as N — > oo. The (normalized) server delivery time is taken as to = max l i=1 {ti} . 
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The normalized number of coded packets user i receives by its delivery time ti over a channel with the 
packet erasure rate e$ is on average U{1 — Let x — 1 — u in ([5]) and oj = ti(l — ej); we have 

(1 - € t )UP\x) + ln(l - a:) > 0, Vx e [0, (6) 

and consequently, 

U = inf{r : (1 - e^rP' (x) + ln(l - x) > 0,Vai e [0,z,)}. (7) 

On the other hand, the recovered fraction by the user on a link of loss rate e as a function of the number t 
of transmitted packets (normalized by N) is 

z(t, e) = sup{C : (1 - e)tP'(x) + ln(l - x) > 0, Va; e [0, £)}, (8) 

i.e., the recoverable fraction is the maximum £ such that the expected ripple size stays positive in the whole 
range of [0, £]. 

We now have all elements to state an optimization problem for minimize the server delivery time t. Using 
(jSJ), the optimization problem is expressed as follows: 

min.p i4l j, t = max U (9) 

i 

s.t. -ei)P'(x) +ln(l -ai) > 0, < x < z h for i = 1, 2, . . . , I; 

or equivalently, 

min.p !to t (10) 
s.t. i (l - £i)P'(x) +ln(l - a;) > 0, < x < zu for i = 1,2,..., Z. 
P(l) = 1. 

Claim 3. There must exist an optimal solution to Problem (|10p tuit/i a polynomial P(x) of degree no higher 
than d max = r i-m^^} ! - I- 

Proof. This claim is provable with an argument similar to Lemma 2 of [14j |. Suppose (t* , P*) is an optimal 
solution of Problem (fit))) , construct P such that pj = p* for j = 1, 2, . . . , d max — 1, and p,i max = 5^,->j x ■ 
Then P(a;) still represents a degree distribution, and meanwhile, 

i'=l J>1 

= E .//';•'•' : • <u. E pj**-- 1 - E : E ^ 1 
= E iM /m " 

j>d max 

= E Pi^ -1 ^-^'-^] 

i>d max + l 

^ E Pj^^Vmax - (dmax + l)^]. 

i>d„ la x + l 

As long as x < ^^T, i.e., d max > J- = ^ - 1, > P*'(x). Thus, set d max = \ - ] - 1 > 

[tzj - ! — 1 > tz 1- ~ 1) we have for i = 1,2, ... ,1, 

t*(l - ei)P'(x) +ln(l- a;) > - £i)P*'(a;) + ln(l - a;) > 0, < x < z h 



G 



and hence (t* , P) is also a feasible and optimal solution of Problem (fTU)) with optimal value t* , and the 
highest degree of P is no more than d max = \ j=^[-^7j ] - 1- 

Thus, Problem (|10[) can readily be converted into a linear programming problem by the method proposed 
14j . For j = 1,2,..., d max , let a, = tpj, and Problem (fTUj) becomes 



mm. 0li ... )(I<w 2^ a 3 (11) 

3=1 

s.t. jajX^ 1 > — — — , < x < Zi, for i = 1, 2, . . . , /. 

i t — £j 

Oj > 0, j = 1, 2, . . . ,d max . 

To solve the linear programming problem (jlip numerically, the constraints defined on a continuous 
interval of parameter x are written out by evaluating x at discrete points within the interval. Lower bounds 
for the minimum value of (TTT|) . and © can thus be obtained. (In Section @J we interestingly observe that 
in a 2-user scenario, the optimal server delivery time obtained from the optimization described here is close 
to the delivery time achievable by using a time-sharing scheme to broadcast degraded message sets. The 
time-sharing scheme is described in 13.4.31 ) 

3.1.4- LT Coding with a Systematic (Uncoded) Phase 

We also study a variation of the LT codes that start with a systematic phase, namely, transmission of all 
the original uncoded content packets (systematic packets) followed by parity packets. 

The formulation of the degree distribution optimization problem is essentially the same, except that © 
(describing the condition which the number of transmissions and the degree distribution must satisfy to 
allow the delivery of the demand of each user i) becomes constraint (fT2|) in the following Claim 0] 

Claim 4. Suppose after the systematic phase, the degree distribution of the parity packets transmitted subse- 
quently follows the distribution represented by P{x), as N — > oo. Then, to recover a fraction Zi, the normalized 
number of parity packets transmitted should satisfy 

-ln(e i ) + (l-e i )(t i -l)P'(x)+ln(l-x)>0, Vx E (1 - e u z,). (12) 

Proof. Please refer to the Appendix. 

The new optimization problem is still readily transformable into a linear programming problem. 
In addition, we have 



1_ e . ■> Zi ^ 1 £i 

inf{V: _ln(e i ) + (l-e i )(r i - l)P'(x) + ln(l - x) >0,Va; £ [0,Zj)}> *i > 1 " e * 



(13) 



The systematic phase delivers the demand of a user with Zi < 1 — which helps reduce the server 
delivery time. It is similarly possible to formulate optimization problems to allow transmitting less or more 
than one round of systematic symbols, and find out the tradeoff between the fraction of systematic symbols 
and non-systematic symbols. However, this is beyond the scope of this paper. 

3.2. Growth Codes 

Growth codes were proposed by Kamra et al. in [7[ to improve data persistence in sensor networks in 
face of sensor node failure. Growth codes were not designed for our heterogeneous scenario described in 
Section[2] However, their feature of progressive partial recovery suggests that they may be a good candidate 



scheme. Extensions and applications of growth codes to video streaming has been studied in, e.g., 16L Il7| 



Particularly, in 16j, a systematic version of growth codes with unequal protection for layer coded video 



content was proposed. 
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Growth coded packets are, as LT coded packets, binary random linear combinations of the original source 
packets, and can be decoded by a belief propagation decoder that is essentially identical to that of the LT 
codes, as described in Section 13.1.11 But, unlike an LT coded stream which produces statistically identical 
packets, a growth coded stream starts with degree-1 coded packets and gradually move on to send coded 
packets of higher degrees. The encoding scheme described in Q operates as follows. Let Rj = —p^j- for 

j = 0, 1, 2, . . . , N - 1. Let Aj = Eif|^_ lJ+ i (°){n- s) for j = 1, 2, . . . , JV - 1. Then, on a perfect erasure 

free channel, the source node sends A\ degree-1 coded packets followed by A 2 degree- 2 coded packets, A 3 
degree-3 packets, and so on. Demand of size between Rj-i and Rj is expected to be fulfilled during the 
phase when degree- j coded packets are sent. 

Growth codes are based on the design philosophy to greedily ensure the highest probability of recovering 
a new content packet upon receiving each additional coded packet. On a link with loss probability e, it is 
reasonable to scale each duration Aj in which degree- j coded packets are transmitted by 1/(1 — e) so as to 
keep the degree distribution of the packets reaching the sink approximately the same as if the code runs on a 
perfect channel. There is no straightforward way to extend such design philosophy and scaling approach to 
broadcasting over channels of different erasure rates, but we can still scale AjS for one of the users and see 
the resulting delivery time for other users and the server, and search for the scaling factor with which the 
server delivery time is minimized. We use the belief propagation LT decoder to decode growth codes, and 
use to predict the delivery time (as defined in Definition [3]) and recoverable fraction of the scaled version 
of growth codes as N — > oo. We compare growth codes with the optimized LT codes in our heterogenous 
scenarios in Section 0] 

3.3. Chunked Codes 

Chunked coding was first proposed and studied in Q. It has also been investigated by the authors of 
this paper in Q via a coupon collection analysis. The TV content packets are grouped into n disjoint subsets 
of h packets (assuming N is a multiple of h), and these subsets are called "chunks". Packets are represented 
as vectors of symbols from GF(g). In each transmission, we first uniformly, randomly select a chunk, and 
then sample from GF h (g) a coding vector uniformly at random and form a linear combination of the content 
packets within the selected chunk. As soon as a sink node has collected h coded packets with linearly 
independent coding vectors generated from the same chunk, all the packets of this chunk can be decoded by 
performing Gaussian elimination on GF(g). As opposed to full network coding, with which coding vectors 
are chosen from GF JV (g), this scheme has lower computational complexity and also to some extent allows 
partial recovery. 

From Theorem 4 of [9j, we know that, for a sufficiently large field size g, the expected number of trans- 
missions needed to decode any k of the n chunks on a unicast link subject to packet loss rate of e is given 
as 

E[T(n,k,e)] = ^ J" |e i^) S t~ 3 ( x ) ^ ~ S h (x)}^ e^dx, (14) 

where 



x 



-l 



^(«)= l +ii + 2T + - + (^i)i (m - 1} 

Soo{x) =exp(a;) and So(x) = 0. 
The above result is based on the generalized birthday problem in [3] ■ 

Definition 4. (Chunked Codes Delivery Time) With the chunked-coding scheme, the normalized delivery 
time of user i is defined as ti = E[T(n, [~-?j7i"|, e^)]/iV. The normalized server delivery time is defined as 
to = max' =1 ti. 
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3-4- Lower Bound and Reference Schemes 

3.4-1- A Lower Bound 

An obvious lower bound for the minimum server delivery time is t — ma,x l i=1 { ^ }. 

3- 4-2. Multiple Unicasts 

If instead of broadcast, the server transmits separate streams to each user, the total normalized number 
of transmissions required will be at least t = ** . 

3.4-3. Braodcast Degraded Message Sets by Timesharing 

Here we describe a reference coding scheme. Without loss of generality, assume z$ = < z\ < zi < 
■■■ _ zi. Then, segment N descriptions/packets into / layers, with Layer i(i — 1,2,..., I) containing 
Li = (zi — Zi—t)N packets. Protect Layer i by an erasure code of rate R4 = 1 — maxjei, e$+i, . . . , e/}, and 
transmit the protected layers sequentially. When N goes to infinity, there exist erasure codes that allow the 
server to deliver Layers 1 through i, that is ZiN packets, to user i for all i = 1, 2, . . . , I in Y^\=i time. 
Later we will find in Section 14.21 that in a 2-user scenario the server delivery time of this scheme is close to 
that of the optimized LT codes without a systematic phase. 

4. Performance Comparison 

In this section we evaluate the schemes described in Section [3] by numerical calculation and simulation 
for a 2-user scenario. 

4- 1- Partial Recovery Curves 

We demonstrate in Figure[2]the evolution of the fraction of recoverable content packets at sink nodes with 
the growth of the number of transmissions from the source in the 2-user broadcast scenario. The erasure 
rates are t\ =0.1 and £2 = 0.5. User 2 has a worse channel. Shown are the performance curves of optimized 
LT codes (with and without a systematic phase), growth codes, and chunked codes at both users. The details 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

t 1 



(a) Numerical Results (b) Simulation Results 

Figure 2: Evolution of recoverable fraction with time at each user in a 2-user scenario, ei = 0.1, €2 = 0.5. N = 1024. Schemes: 
(1) LT optimized with z\ = 15/16 and Z2 = 9/16. (2)optimal systematic LT (3) growth codes with a scaling factor of 1 _ L e 
(4) chunked codes with n = 16 chunks of 64 packets. LT(-sys): Optimized LT codes (with a systematic phase); GC1: Growth 
codes scaled by j-jj- ; CC64: Chunked codes with 64 packets per chunk. 

in obtaining the numerical results are listed below. 
• LT codes (systematic and non-systematic) 
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— The degree distributions are obtained by solving the optimization problems in 13.1.31 and 13.1.41 
by setting (zi f ei) — (15/16,0.1) and (22,^2) = (9/16,0.5). For the non-systematic version, the 
optimal P[x) = 0.0195a; + 0.7814a; 2 + 0.1991a; 3 ; for the systematic version, the optimal P{x) — 
0.7061a; 2 + 0.2939a; 3 . 

— The delivery time is computed for given Z\ and based on ([7]) in the nonsystematic case and 
on ([T3"l) in the systematic case. The optimal server delivery time is 1.5178 for the nonsystematic 
version, and 1.2488 in the systematic version. 

• Growth codes 

— N — 1024. The time spent transmitting coded packets of each degree i, Ai, is scaled by a factor 
of 1/(1 — ei), that is, the code is adapted to the channel conditions of user 1. 

— The recoverable fraction is computed by evaluating (JSJ on the empirical degree distribution of 
received packets. Note that (JSJ) assumes N — >■ 00. 

• Chunked Codes 

— We set the number of packets to N = 1024, the number of chunks to n = 16, and thus the chunk 
size becomes h = N/n = 64. 

— We compute E[T(n, \zn\,€i)]/N with E[T(n, k, e)] given by (fT4| . 

Simulation results show the average time required for the user to recover a given fraction z of all the 
packets. The average is taken over 100 runs, and the number of source packets N is set to be 1024 in all 
runs. For chunked codes, the finite field size q is set to be 256. It is shown that the simulation results are 
close to the theoretical performance prediction for the coding schemes studied, except for the difference in 
the latter stage of growth codes, which comes from the slight difference in the coding scheme numerically 
evaluated and the one simulated. In our simulations, after finishing sending A m degree-m, coded packets 
where m = [~ max . ^{l-z-} ! > instead of proceeding to sending packets of higher degrees, as defined in the 

original scheme, we allow the server to send coded packets according to a degree distribution pj = ^„"f 3 , 
for j — 1, 2, . . . , m, the cumulative degree distribution up to stage A m . This provides the user on the worse 
channel with the low-degree packets that are needed but have been lost due to a higher packet loss rate, and 
allows the user to proceed in packet recovery. All the coding schemes exhibit partial recovery property to a 
corresponding degree, namely, the recoverable fraction gradually increases with time. With both optimized 
LT (without the systematic phase) and the chunked codes, there is an initial stage where no packets are 
recoverable. 

4-2. Delivery Time 

We continue with a study of the server delivery time performance of the coding schemes in the 2-user 
scenario. We keep e\ — 0.1 and £2 = 0.5, set zi — 15/16, let Z2 vary from to 15/16, and plot the server 
delivery time versus z 2 . With LT codes, for each set of (zi,z 2 ) values, we solve for the optimal degree 
distribution. With growth codes, for each Z2, we search all scaling factors between [ 1 _ e , for the one 

that minimizes the server delivery time. With chunked codes, we examine all power-of-2 chunk sizes between 
h = 2 = 1 and h = 2 10 = iV, and find the chunk size that minimizes the server delivery time. The numerical 
evaluation of the minimized server delivery time is plotted versus the demand Z2 of user 2 in Figure G3 Note 
that the optimal degree distribution changes as the demands vary. 

The server delivery time of the reference scheme via time-shared broadcast of degraded message sets, as 
described in 13.4.31 is also plotted. In the two- user case discussed here, z\ > z-i and t\ > €2- Therefore, two 
layers are formed, the first consisting of L\ ~ z 2 N blocks and the second another L 2 — (zi~z 2 )N blocks, and 
respectively protected with erasure codes of rates 1 — e\ and 1 — e 2 ■ It is particularly interesting to find that in 
the setting demonstrated here, the reference scheme performance almost coincides with the non-systematic 
optimized LT coded scheme with optimized degree distribution. The reason of such a phenomenon, however, 
awaits further investigation and is beyond the scope of this work. 
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Figure 3: Comparison of the server delivery time t: N = 1024 descriptions/packets (N — > oo for LT codes and growth codes), 
ei = 0.1, €2 = 0.5, 2i = 15/16, 22 ranging from to 15/16. 



In addition, in Figurc[3l we plot two other reference lines: t — max! . Zl . -r^ 2 — j and t — -r^ 1 — h , Z2 . The 

' ° 1 — * " L 1 — ei ' 1 — 62 J 1 — ei 1~£2 

former represents an apparent lower bound for server delivery time, and the latter represents a lower bound 
for the total number of packet transmissions if the demands are delivered by unicasts to each individual user. 

Comparing the server delivery time of systematic and nonsystematic LT codes optimized with the knowl- 
edge of channel state and user demands, we find that the systematic phase significantly shortens the server 
delivery time when the demand of the user on the worse channel is low. As the demand of user 2 rises, 
however, including a systematic phase becomes suboptimal, until the demand rises further closer to 1, when 
the delivery time seems to converge to the delivery time with the nonsystematic scheme. 

Growth codes and chunked codes do not have a competitive server delivery time performance compared 
with other schemes. When the demand of the user on the worse channel is low, using these schemes requires 
more transmissions to deliver the demands than using multiple unicasts. One advantage of chunked codes, 
however, is that all the content packets in the same chunk are decoded simultaneously, allowing these packets 
to be dependent on each other for the reconstruction purpose, which may reduce redundancy in the source 
coding stage. 

It is also of interest to study performance measures other than the delivery time of the server. For example, 
one could wish to maximize the minimum of users' throughputs rninj{zi/ti}, the ratio of the demand and the 
time needed to complete the demand, or to maximize the minimum channel utilization mini{^/((l — e^)^)}, 
the ratio of the information pushed through the channel to the channel bandwidth. These criteria are of 
interest for acieving, e.g., fairness in the network, and optimization for these criteria can yield very different 
code design parameters (degree distributions). The optimization of the degree distributions for LT coded 
broadcast for these different criteria has been studied in our work f5~oj ] and interested readers are kindly 
referred to this paper for detailed problem formulation and results. An additional type of heterogeneity 
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can also be treated in the framework provided in Section 13.11 namely, when some of the sink nodes cannot 
decode but can only recover content from degree-1 coded packets. Readers are also referred to [l(| for results 
regarding the coexistence of nodes able and unable to decode in the system. 

5. Conclusion and Future Work 

We investigated the usage of three coded schemes for broadcasting multiple description coded content 
in a single-hop wireless network with heterogeneous user nodes of nonuniform demand over links of diverse 
packet loss rates. The three coded schemes include the LT codes with specially optimized degree distribution 
(with or without a systematic phase), growth codes, and chunked codes. Particularly, with the LT codes, 
we are able to formulate the degree distribution design problem in the heterogeneous scenario into linear 
optimization problems. For the schemes compared, we characterize numerically the number of transmissions 
needed to deliver the demand of all users. The numerical evaluations agree with simulation results. 

A systematic phase delivers the demands efficiently when the fraction of content requested by the users 
does not exceed the link capacity. On the other hand, for higher demands, coding significantly improves the 
delivery time. Different user demographics result in very different coding schemes being suitable for efficient 
delivery of demands. Growth codes and chunked codes, are not found to be as suitable to the communication 
scenario as the optimized LT codes. Interestingly, time-shared broadcast of degraded message sets is found 
to give a comparable delivery time performance as that of the non-systematic optimized LT codes. 

As for future work, we are interested in incorporating the source coding stage into code design, since 
there is clearly an interplay between the compression efficiency in source coding and the efficiency of channel 
coding. We are also interested in the exploration and analysis of more competitive schemes for broadcasting 
to heterogeneous users. 
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Proof of Claim [2] 

We give two independent arguments to validate the claim. 

After the systematic phase, each content packet successfully reaches user i independently at probability 
1 — ej, and an average of (1 — ei)N packets reaches the user. If Zi < 1 — e,, the demand of user i is considered 
fulfilled. We consider the case where Zj > 1 — ej. 

In the first argument, we find out the non-systematic equivalent degree distribution P(x). To do that 
we only need to find out the (normalized) number of degree-1 packets(containing content packets selected 
uniformly at random with replacement) required to be transmitted in order for user i to recover the (1 — 
€i)N distinct content packets received in the systematic phase. This number, — , can be obtained by 
setting P(x) = 1 • x(all-one degree distribution) in (J7J. A coupon collector's argument brings to the same 
conclusion [19], Ch. 2] (see also [20(): the expected number of samplings required to collect zN distinct 
coupons is 



V 



/ll 1 W N ( zN- 1\ 

( — + + ••• + ) > Win = -Vlnf 1 J. (.1) 

\N N-l N-zN+lJ- N-zN+l V N J y ' 



Divide (TT|) by N(l — e,), take V — > oo and let z = 1 — ei, we find the same result as found from ([7]). 

Hence, pi = -£rj + Pi and P(x) = -jz^x + P(x). In ©, replace P(x) by P(x) and U by (U - 1) to 
subtract the time spent in the systematic phase, and we get (TlT 
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Alternatively, one can consider the packets received in the systematic phase as side information to the 
decoder, such as in [? ? ]. These packets are removed from the subsequent coded packets, and we need to 
decode another (z, — (1 — 6i))N packets from the remaining e^iV packets. Let d and d be the random variables 
representing the degree of a randomly generated coded packet and the degree of the coded packet after the 
removal of the packets received in the systematic phase. Thus, d = Ylj=i Xj where Xj(j — 1, 2, . . . , D) arc 
i.i.d Bernoulli(ei) random variables indicating if the jth content packet participating in the coded packet has 
not been received in the systematic phase and remains in the coded packet. Thus, the moment generating 
function D is P{x) — P(Q(x)) — P(l — t{ + eix) where Q(x) = 1 — ej + tiX is the moment generating 



function of the i.i.d. 



X jS . 



In ([B]), replace P(x) by P(x), Zi by 



t, 



_ (tj-l)N _ ti 



( Zl -(l~£,))jV 



we have 



-(i-^) 



and ti by 



U - 1 



(1 - ei)- P'(Q{x))Q'(x) + ln(l - a;) > 0, Vx e 



0, 



Zi - (1 - £,*) 



(■2) 



Let y = Q(x) = 1 — 6j + 6iX, we get 



(1 - ,;u'—^l", !r . + In (l 11 ' ' '' ] v o. v,, :: 



1 - e,-. 1 - 6i + e, 



^ - (1 - Cf) 



(.3) 



which is exactly (|T2"j) . 
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